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Abstract 

We present exact solutions in Einstein- Yang-Mills-Dirac theories with 
gauge groups SU{2) and 5[/(4) in Robertson- Walker space-time IR, x 5^, 
which are symmetric under the action of the group 50(4) of spatial ro- 
tations. Our approach is based on the dimensional reduction method for 
gauge and gravitational fields and relates symmetric solutions in EYMD 
theory to certain solutions of an effective dynamical system. 

We interpret our solutions as cosmological solutions with an oscillating 
Yang-Mills field passing between topologically distinct vacua. The explicit 
form of the solution for spinor field shows that its energy changes the sign 
during the evolution of the Yang- Mills field from one vacuum to the other, 
which can be considered as production or annihilation of fermions. 

Among the obtained solutions there is also a static sphaleron-like so- 
lution, which is a cosmological analogue of the first Bartnik-McKinnon 
solution in the presence of fermions. 



1 Introduction 

Exact solutions in gravity coupled to fields of different types have always attracted 
much attention. In particular, the last few years witnessed a great interest in 
solutions to EYM systems. There were found both numerical |]1], 0] and exact 
solutions with 5*0(3) and 5*0(4) groups of spacial symmetry 0, The exact 
5'0(3)-symmetric solutions turn out to be static and singular and are, in fact. 
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a generalization of Reissner-Nordstrom solutions to non-abelian gauge theories. 
The S'0(4)-symmetric solutions correspond to the Robert son- Walker ansatz for 
the metric and are interpreted either as wormhole solutions in the euclidean 
domain ||^ or as cosmological solutions for the radiation dominated universe 0. 

There were attempts to accommodate spinor fields in this picture as well, but 
they usually amounted to considering spinor fields in the EYM background. For 
example, in 0, |], ^ an anomalous fermion production in the EYM background 
was discussed. 

The study of exact solutions in Einstein-Dirac systems also attracted much 
attention [0, |ll|. It was found that the Robertson- Walker ansatz in such a 
system leads to the so-called "ghost solutions" in ED systems, for which the 
Dirac field has a vanishing energy-momentum tensor. 

In the present paper we consider a self-consistent EYMD system and find exact 
solutions for the case of the gauge groups SU{2) and S't/(4). The group SU{4) 
is the simplest gauge group which gives qualitatively new results in comparison 
with the group SU{2). Our solutions describe the radiation dominated universe, 
in which an exchange of energy between the YM and Dirac fields takes place. 
They can be of interest for treating the dynamics of the early universe, because 
they are based on the energy-momentum tensor derived from the fundamental 
Lagrangians of particle physics, rather than on the phenomenological one. Our 
solutions also solve the problem of the "ghost solutions" in ED systems. 

To find these solutions we employ the dimensional reduction method for grav- 
itational, gauge and spinor fields |T2|, [T^, which enables us to relate symmetric 
solutions in EYMD system to certain solutions of an effective dynamical system 
and essentially simplifies the problem of finding symmetric solutions. 



2 The effective action 

We consider an 5'0(4)-symmetric Einstein- Yang-Mills-Dirac system with the 
standard action in space-time M = R x S''^: 

S = Se + Sym + Sn , (1) 

where 

Se = [ (R-A)dz;, (2) 

IQttk Jm 

Sym = ^ [ tr(F,^r^)dt;, (3) 

Sd = ^J^7"i?a'^(^-^^M-/5[7^7^]+m^))v^d^ + /^•c.. (4) 

Here A denotes the cosmological constant, R is the scalar curvature and dv = 
■\/|det(g)|d'^x is the canonical volume form corresponding to a metric g on M = 
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R X S"^. F and A are the Yang-Mills field strength and potential respectively, 
and the trace is taken in the adjoint representation of the Lie algebra of the 
gauge group. In the Dirac action, denotes the spinor field, {Ea = E^^d^} is 
an orthonormal frame on M = R x 5^, {7"} arc the gamma matrices and 5 is 
the representation of the gauge group in spinor space. 

We identify with the group manifold of SU{2). Then the action a of 
K = ,50(4) = {SU{2) X SU{2))/{1, -1} on SU{2) is given by 

a{{h,k2),x) =hxk2~\ {h,k2)eSU{2)xSU{2),xeSU{2). (5) 

The isotropy group H at x — Isu{2) is isomorphic to SU{2) and is given by 

H = {{k,k)eK, keSU{2)}. (6) 

Further, we have the reductive decomposition ^ — (B^Xft of the Lie algebra ^ 
of K, where 

Sj = {{X,X); Xesu{2)}, (7) 
m = {{X, -X) ; X e su{2)} . (8) 

Obviously = {SU{2) x SU{2))/H is a symmetric space, and S'0(4)-invariance 
is equivalent to invariance under SU{2) x SU{2). The isotropy representation of 
H — SU{2) in 271 is 3, i.e. it is isomorphic to the adjoint representation. 

Now we have to reduce the action of the original theory due to the 5*0 (4)- 
symmetry. We begin with the action of the gravitational field. 

The most general form of an 5'0(4)-invariant metric g on M is 

g = -N^dt^ + a^dnla , (9) 

where dfl^ga is the standard metric on a 3-sphere of radius 2 resp. — | the Killing 
metric on SU{2) = S^. For the sake of the future convenience we assume N and 
a to have the dimension of length and t and dl^la to be dimensionless. 
An orthonormal coframe {0'^} on M is given by 

e° = Ndt, 6' = a-&\ i = l,2,3, (10) 

where '(9Vj/(2i) = i) is the canonical left invariant 1-form on SU{2), ai being 
the Pauli matrices. In what follows we set = ak/{2i) and denote by rjap the 
Minkowski metric rj = diag(— 1, 1, 1, 1). It is a matter of simple calculations to 
find the components of the spin connection on M: 

i^ik = -^£ikiO\ i=l,2,3. (11) 
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Using the standard formulae for the curvature, substituting it into (0) and inte- 
grating over S^, one easily finds the reduced gravitational action 



where IGvr^ is the volume of SU{2) = with the standard metric di7|3. Omitting 
a complete divergence, we get the effective action 

^1^ = 1?- f (^aN-Aa'N-6^aa')dt, (13) 

which we consider as the reduced action of the gravitational field. 

Next we turn to the gauge field action (^. An 5'0(4)-symmetric gauge po- 
tential A on R X S*^ is in one-to-one correspondence with a triplet {t,A,^}, 
where r is a homomorphism from the isotropy group H into the gauge group G 

t:H^G, (14) 

A is a gauge potential on R with values in the centrahzer ft(^{T'{S))) of t'{S^) in 
= Lie{G) and $ is a linear mapping 

$ : R ^ 971* ® (5 (15) 

with 



$oAd(/i) = Ad(r(/?.))o<l>, \fh e H . (16) 

Here Ad{h) is the restriction to H of the adjoint representation of K applied 
to OK and Ad(r(/i)) is the restriction to t{H) of the adjoint representation of G 
in (5. As we have already mentioned, we will consider gauge groups SU{2) and 
SU{A). In the first case the centralizer is trivial, and there is no reduced gauge 
potential A, whereas in the second case there can be a nontrivial centralizer. 
Therefore, the case of the group SU (4) is more general, and we will carry out 
all the calculations for this case and then explain the difference from the case of 
SU{2). 

Thus, we take the gauge group G = SU{4) and define the homomorphism 
T : H ^ G hj the decomposition of the fundamental representation 4 of SU(4) : 

4^(2,2), (17) 

i.e. we represent the space C"^ by the tensor product ® and let the fun- 
damental representation of t{H) act on the first factor and the fundamental 
representation of the centralizer Gg{j{H)) = SU{2) act on the second factor. 
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It is known that the adjoint representation of sl{n) can be expressed in terms 
of its fundamental representation by 

ads/(n) = n*^n, (18) 

where tilde means dropping out a one- dimensional trivial representation and n* 
denotes the contragradient representation, t{x)* = —t{x)'^. Therefore we get 

adsM(4) = (2, 2)*® (2, 2) ^ (3, 1) © (3,3) © (1,3) . (19) 

A basis in sm(4) adapted to this decomposition is 



Hk 


= Tk 


g>I = 


PkA 


= Tk 




Ha 


= Ig 


) Ta , 



Tk 
Tk 



(20) 

(21) 
(22) 



where equation ( pof ) explains our rule for evaluating the tensor product of matri- 
ces. It is easy to calculate the commutators of these generators, for instance 

[PkA, Pib] = £klmHm5AB + ^ABcHcSkl ■ (23) 

The constraint ( [T6| ) means that the mapping $ is an intertwining operator, which 
intertwines the isotropy representation 3 of if = SU{2) in 971 with the represen- 
tation (p!9|). We introduce basic intertwining operators from 971 into (& defined 
by the relations 

n{rk,-rk)) = Hk, (24) 
lAiirk,-Tk)) = PkA. (25) 

Then for any X = (X^Tk) G su{2) we have 

$((X,-X)) = ^IiX,-X)+^^lA{X,-X) (26) 
= X''{^Hk + ^^PkA) = X©(P + eVA), (27) 

where ^ and ^^,A = 1, 2, 3 are real valued functions on R. In what follows we 
denote the vector (^\^^,^^) by ^and set ^ = P ^^a^ and | = ^^cr^i. The 
centralizer Cs„(4)(r'(sM(2))) of r'(su(2)) in (S = su(4) is su{2) and is spanned by 
the Lie algebra elements Ha, A = 1,2, 3. The matrix is in the representation 
(1 + 3) of the centralizer (t and e ^^(2). 

The symmetric gauge potential A on R x S"^ can be easily expressed in terms 
of the matrix ^ and the canonical left invariant 1-form ■(9 = ^9*rj on SU{2) = S^: 

A = ^H,^' + ^^^PkA^' + A^HAdt 

= ® + l) + I © Aodt . (28) 
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Here Aq is a function on R with values in su{2), i.e. Aq = AqTa- The term 
I ® Aodt is the reduced gauge potential on R, which can be gauged out, but 
we keep it for the moment, because it is necessary for deriving the equations of 
motion. It is not difficult to calculate the corresponding field strength F and to 
obtain the reduced Yang-Mills action, which takes the simple form 



Sym = — ^ / 24 f H (f+AnxO 



8a 



{{e+{i?-i?+^e{^r) ]dt, (29) 



where Aq is the vector (Aq, Aq, Aq) and Aq x ^ is the vector product of Aq and ^. 

If the gauge group is G = SU{2), the unique nontrivial homomorphism r : 
H ^ G can be defined by the identity mapping, i.e. 

T{{k,k)) = k, keSU{2). (30) 

In this case the centralizer C(g(r'(i3)) is trivial, and the intertwining operator is 

$((X,-X))=e(t)X, XGsn(2),e(t)GR. (31) 

It is clear that the gauge potential A on R x S*^ still has the form (^) with 
,^ = 0, Aq = 0. The reduced Yang-Mills action is also given by (PD|), provided one 
puts ^ = and rescales the coupling constant g 2g. 

Next we have to reduce the action (^) for symmetric spinor field. We choose 
5 to be the fundamental representation, i.e. we can write the spinor '0 as a 4 x 4 
resp. 4x2 matrix on which an element g of S'f/(4) resp. SU{2) acts via right 
multiphcation by g~^. 

In accordance with our choice of the metric signature (H), we have 

{Y,Y} = -2r\ V'" = dzagi- + ++) , 7^ = Z7V7V , (32) 
.0 _ / I \ f -a' \ . / I 



^ =U 0j'^^=U^ J'^ =l0 _iJ^-M,3. (33) 
The bispinor representation A(s) is defined by 

A(s)-S'^A(.) = A(s)^^7^ (34) 

where s is an element of the group Spin{l, 3) and A is the covering homomorphism 
from Spin{l,3) onto 5*0(1,3). On Lie algebra level we get 

[A'{A),Y] = -A^.Y , A'(A) = -^A^[7^, 7l , (35) 
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where A^,^ = —Ay^ is an element of so(l, 3) = spin{l, 3) = sl{2, C). 

An S'0(4)-symmetric spinor field on R x S''^ is in one-to-one correspondence 
with a matrix valued function p on R, which satisfies the condition 

{6'{T\h)) + A'{\\h)))p = ,WheS). (36) 

Here X' : Sj ^ so{l, 3) is the homomorphism induced by the isotropy representa- 
tion, which can be calculated explicitly: A'(ra)^^ = —Sac- Therefore, if the gauge 
group G is SU (4) equation ( ^6]) reads 



I'l'e,,^ + 5' (iJfc) \p = {{r^®\)p-p (r, ® I)) = . (37) 



The general solution of this constraint equation is 



fll / \ Vi V2 

i.e. a symmetric spinor on R x S''^ is parameterized by two complex doublets 
u = {ui,U2)^ and v = (fi,t>2)^, one for each chirality. We see from (|38D that 
the reduced gauge group C = SU{2) acts on both doublets by the fundamental 
representation. Taking into account equations (0) and (p8|) , it is a matter of 
simple calculations to get the reduced action 



Sd = IGtt^ / a^N ( —{u u — uu + vv — vv) 



1 3 " " \ 

+—A^{uaBU + vaBv) - —{u^u-v^v)jdt. (39) 

If the gauge group G is SU{2), equation (^) reads 

{{Tk(^I)p-pTk) = 0, (40) 

and we get 

P=(^J)=I®(^) > u,veC-{R), (41) 

i.e. a symmetric spinor on R x 5*^ for G = SU (2) depends on two arbitrary 
complex functions u and v, one for each chirality. The reduced action has the 
same form (^9]), if we put there ,^ = 0, Aq = 0. 

Now we can write down the reduced action of the coupled EYMD system. In 
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the case of the gauge group SU (4) it has the form 
S — Se + M + Sd 

^ IGtt^ / l^f^aN-Aa'N-Q-a) 

+ ^ia^('U ii — uu + vv — vv) 

+a^AQ{uaBU + vasv) - ^^^{uiu - viv)^ jdi . (42) 

If we choose the gauge group to be SU{2), we have to put ^ = 0, Aq = in this 
action, to rescale the couphng constant g ^ 2g, and to take into account that 
the variables u and v are no longer isospinors, but ordinary functions. 



3 The field equations and solutions 

— * 

Variation of this action with respect toa,N,^,^,u,v and u , v gives us the 
field equations. When taken in the special gauge Aq = 0, a = (the latter 
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condition means that t is now the conformal time), they have the form 
a-^S = ( -a^ - SAa'^ - 6a + 12ad] 

+3a^ (^i{u u — uu + vv — vv) — {u^u — viv)^ = , (43) 



a-^S = -^(h^-Aa^ + Qa' 



5N IQ-KK \2 

3 /:2 1 



[i + {if + ^ ((e^ + {if - If + 4e^(o^ 



2g^ 

-la\uiu-viv) = Q, (44) 



^5 = -^,{2i+{i{e + {if-i)+mf 



Si 2g 



3 

--a^{uu -vv) = 0, (45) 



■^a^{ua^u-va\) = 0, (46) 



(47) 



1 

2^" 


"^^5 


d 3 
= i—a'^u — 
dt 


3.' 3 

-ta^u = 
4 


0, 


1 

-a 
2 


"^^^ 

ou 


. d 3 

= i—a'2u + 
dt 


3 3 

-a^ut = 
4 


0, 


1 

2^" 


6v 


d 3 
= i—a'^v + 
dt 


3.^ 3 

-ta^v = 
4 


0, 


1 

2^^ 


6v 


d 3 
= i—a'^v — 
dt 


3 3^.^ 

-a^vf = 
4 


0. 



Variation with respect to Aq gives a constraint 
6S 3 



(4J 



(sBCDfi^) + a^{uaBU + vaBv) = , (49) 



which means that the total isospin of the gauge and the spinor fields equals zero. 

Now we will show that it is possible to find exact solutions to this system of 
equations. We begin with the simpler case of the gauge group SU{2). There is 
no constraint (^) in this case, and we also have to drop the equation (^Bj), to 
put ^ = in the others and to rescale g 2g. 

The Dirac equations (E^) and (^8]) give 



uu = — , vv = — (50) 
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and 



3i 

uu — uu = ——^uu, (51) 
3i 

V V — V V = —^vv , (52) 

where Cu and C„ are arbitrary positive constants. We will discuss the meaning of 
these constants later, here we assume that they are proportional to the number 
of fermions with positive resp. negative chirality on 5*^. 

Now equations (H3) and (03) simplify considerably. Their sum gives 



^— (Sa^ - 4Aa^ + 12aa) = . (53) 
Multiplying this equation by a/a and integrating we see that 

^a' - Aa^ + 6a' = E , (54) 

where E is an arbitrary constant which has the meaning of the total energy of 
the system. Equation (|5^ is the standard Friedmann equation for the radiation 
dominated universe and has a simple analog in mechanics. We can consider a 
as the coordinate of a particle with mass 1 and energy E/12 which moves in a 
potential 

W^(a) = \o' - ^Aa' • (55) 

If A < Ae = then the motion will be periodical. This means that our 
solution describes a universe which first expands and then contracts, where a = 
corresponds to a singular metric in the beginning and the end. If A > Ae, then 
the solution can be either static or can describe an expanding universe. 

Next we turn to the YM equation (|^) (we recall that we have rescaled g 2g 
in the case under consideration). Due to equation (|3), it decouples from the 
equations for u and v: 

^ [2I + [ae - 1))) + ^(a - a) = . (56) 

The first integral of this equation is 

^(eVi(^-l)^)45(C.-C„, ^ ^, (57) 



where the integration constant is due to equations (|4^) and (| 
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Equation ( pT] ) also has an analogue in mechanics. A point particle with mass 
1, energy g'^E / {127T k) and coordinate ^ moves in a double-well potential 

ViO = l{e-ir + 2g'aC^-C^). (58) 

We can interpret the first term in (|57|) as the energy of the Yang-Mills field and 
the second term as the energy of the Dirac field due to the interaction with the 
gauge field. The equation describes an exchange of energy of the two fields, 
the coupling between the gauge and the spinor field being proportional to the 
difference of the numbers of left and right handed fermions — Cy. The exact 
solution of equation (E% is possible in terms of elliptic functions of the first kind 



We consider for instance the case, where the energy of the system is 
larger then the local maximum of the potential V{C,). Then the system will move 
between the turning points defined by the real zeros of the polynomial 

^-nO- (59) 

In the case under consideration we have two real zeros a and P, a > > P, and 
two complex conjugated zeros m — in and m + in, i.e. 

lia - 0(e - - + = ^ - = \i ■ (60) 
We have to solve the integral 

m = r ■ (61) 

|£-i(^^-i)^-Va:(a-a) 

Its exact solution is given by 



2 -LrccteJ^.i. -<^-^'^ + <"-'^>^ 



. _P |2arcctgy^.-y 1 . (62) 



where F is the elliptic integral of the first kind and p, q are defined by 

= {m- af + , = {m - pf + n'^ . (63) 

This function can be inverted, and we can get ^(t) expressed in terms of the 
Jacobi elliptic function sinus amplitudinis. 

With a given solution ^ we can solve the Dirac equations (^) and ( ^8]) 



u = y^exp|-z jl^dt}, (64) 
^"expjz [^^dt\. (65) 




a3 " / 4 
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One easily checks that our solution fulfills the whole system of field equations. 

Now we consider the case of gauge group SU{4:). The Dirac equations (^) 
and (^8]) again lead to (0), but instead of (^Tj) and (|52|) we now have 

uu — uu = ——uC,u, (66) 

vv — vv = ^viv . (67) 

These equations are also sufficient to decouple the Friedmann equation from the 
Yang-Mills-Dirac equations, and we get again the equation (0) for the scale 
factor a. 



Now we have to solve the Yang- Mills equations ( ^5]) and (^61). We start with 
the discussion of the constraint (|49|). In what follows we restrict ourselves to the 
case, where 

eBCDf^'' = for 5 = 1,2,3, (68) 

i.e. the isospins of the gauge field and the Dirac field vanish separately. This 
equation means that the angular momentum of the motion in the ^-space equals 
zero, that is the motion goes along a straight line passing through the origin, and 
the vector ^ is always proportional to a fixed vector ^o- We use the remaining 
gauge freedom to choose = (0, 0, 1) , i.e. ^ = (a^,, C € R. In this gauge we 
obtain from the Yang-Mills equations (^) 

u(TAU — vaAV = for A = 1,2. (69) 



Further we get from equations (^9|) and (|6q) 

uctbu + vctbv = , B = 1,2,3. (70) 

Hence, we have uaiu = ua2U = 0. These equations have two solutions: 

1st case: u = a~^^'^aiw^ ,v = a~^''^/5iw_ (71) 
2nd case: u = a~^^'^a2W- , v = a~^^'^(32W+ , (72) 

where w+ = (^^^^ , w_ = (^^^^ (73) 

and ai, Pi e C, \ai\ = = 1,2. In particular, we obtain Cu = = |aip 

in the i'th case. If ^ = ^^"3 = 0, then the Yang-Mills equations (^6]) demand 
u = V = 0. Therefore, we consider only the nontrivial case, when ^ ^ 0. Hence 
we have two Yang-Mills equations (^) and (^6]), which now take the form: 

= o, (74) 
^{2c+{c{e+e-i)+2ec))+^s = o, (75) 
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with S := Cu 



\ai\ 



in the first and S := —Cu = — |«2p in the second case. 



To solve equations ([7^) and (^), we pass to new variables in accordance with 



(76) 



It is easy to check that the equations for x and y decouple and take the form 



X[X 



1))+3S 



The first integrals of these equations are 



3 
2? 



. 2 



1 

4' 



1)^ 



3S 



3Sx 



3Sy 



0, 
0. 

El 
E2 



where due to 



and (0) the constants E fulfill 
El + E2 = E . 



(77) 
(78) 

(79) 
(80) 

(81) 



These equations can also be solved exactly in terms of Jacobi elliptic functions, 
but unlike the case of the group SU{2), there will be two different periods of 
motion in x and y. 

We would like to note here that equations of the type (0), ( [75D with S = 
for the Euclidean EYM system were first found in ||^ , but they were solved there 
only for the case either ^ = or = 0. 



Substituting the solutions for ^ and ( into equation (p8|) we get 



A 



^^e 

l+VQ 







^2) 



i.e. the gauge potential A takes values only in an su{2) © su{2) subalgebra of 
sm(4), each su{2) part of the gauge potential being coupled to only one of the 
spinor fields u resp. v. 

If we have a solution to the system of equations (|79D and (|80D it is easy to 



integrate the Dirac equations 
obtain 



u 



and 



— exp 

Cu 

exp 



With given solutions x and y we 




3 , 
- xdt 
4 



ydt 



^3) 
M) 
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in the first and 




in the second case. 



4 Discussion 

In studying the self-consistent EYMD system we found that the evolution of the 
metric decouples from the remaining system and is described by the Friedmann 
equation for the radiation dominated universe. The same result for the case of 
EYM systems was obtained earlier in ^, and solutions for the spinor field in 
this background were studied in 
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Unlike the latter solutions, our solutions 
take into account the back reaction of the spinor field on the YM field. 

The Yang- Mills equations in the case G = SU (2) admit three static solutions: 
two minima and one local maximum of the potential 



(87) 



Of course, these are solutions of the whole system only if the constant E is equal 
to V{^). If we have Cu = C^, we can interpret these extrema as two vacua 
(^ = — 1,^ = +1) with Chern-Simons numbers and 1 and as a sphaleron-like 
solution = 0) with Chern-Simons number 1/2 lying on top of the potential 
barrier between the vacua 0, ^. Chiral spinor fields shift slightly the location 



of the extrema and the Chern-Simons index of the corresponding gauge field 
configurations, which is quite natural in the presence of matter fields [O, |18|, f9 



By fine tuning the cosmological constant A and the energy E, we can get a 
static sphaleron-like solution of the whole system. This solution corresponds to 
the local maxima of V{^) and W{a), see equations ( p5D and Obviously, this 
solution has two unstable modes - one in the gravitational and one in the gauge 
field sector. This is another indication that the static solution is a cosmological 



analog of the first Bartnik-McKinnon solution |20, 21, 22 



At this point, we have to comment on the meaning of our classical spinor 
field. The problem of interpreting spinor fields in cosmology has been discussed 



for many years pOl , |rT[, but so far no satisfactory solution has been found. It is 
clear that our solution describes just one energy level of the Dirac field in the 
EYM background, which is exactly the so called zero mode. If we assume that 
the Dirac field is normalized to unity, the influence of this field on the EYM 
system is negligible. Therefore, we suggest that the Dirac field of our solution 
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is normalized arbitrarily and is, in fact, an effective field describing fermionic 
matter with S'0(4)-invariant energy- momentum tensor. This assumptions seems 
to be reasonable in the cosmological setting, because the energy levels of the Dirac 
field in the EYM background must be very dense, and replacing the contribution 
of fermions on the lowest levels by that of the zero mode level could be a good 
approximation. 

Equation (|4^ ) is the (0, 0) component of the Einstein equations, and therefore 



we can interpret the constant E as the total energy of the YMD system. On the 
other hand, equations ( ^7|) resp. ( [79l) and (|80D describe the exchange of energy 



between Yang-Mills and spinor field. This interpretation is supported by the 
solutions (1^) and (^) for the spinor field in the case of the gauge group SU{2): 
the momentary frequency resp. energy of the spinor field is given by the integrand 
in the exponent of the solutions and this is up to a factor ^. 

This observation also means that our solutions describe creation and annihi- 
lation of fermions. If the total energy of the YMD system is larger than the local 
maximum of the potential V{^), the motion in the variable ^, stemming from 
the YM field, will be periodical. When ^ crosses the maximum of the potential 
V, it changes its sign. With the above interpretation we see that the energy of 
the spinor field also changes its sign. This is in accordance with the observations 
in P, P, where it was shown that the spinor field has zero modes in the 
sphaleron-background and that moving between neighbouring vacua of the gauge 
field results in a shift of the energy level of the spinor field. We also obtain from 
our solution that the effect is opposite for the spinor field with opposite chirality. 
The corresponding violation of the fermion number can be calculated explicitly. 

The classical U{1) vector current of the Dirac field is given by 



J = ^7«^^" = ^^^e' . (88) 



This current is classically conserved, i.e. d{*j) = 0. Here the star denotes the 
Hodge star and {6"^} is the orthonormal coframe on R x S^, see equation (|10D . 
But if one considers chiral matter on quantum level, this current has an anomaly 



24, 19 : 



d{*j) = (l/167r2)tr(F A F) . (89) 
Integrating (|89| ) over I x S^, I = [ti,tf], we get 

NF{tj)-Np{U) = {l/167r') [ tr(FAF)=/ Q{tf) - [ Q{U) , (90) 

JlxS^ Js^ 



where Q is the Chern-Simons 3-form 



Q = T7^tr(A A rfA + A [A, A]) , (91) 
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and 

Npit) = [ *J (92) 

J{t}xS3 

is the total fermion number at time t. Hence, for chiral fermions one obtains on 
quantum level that the change in the fermion number is given by the difference 
of the topological charges of the gauge field at the two ends of the time interval 
/. The topological charge 



q= / Q 

at any fixed t can be calculated from the formula for the symmetric potential A, 
which gives 

If we pass from the configuration with ^ = — 1 at t = tj to the configuration with 
^ = latt = t/, the topological charge changes by 1. Setting in our solution one of 
the constants, Cu or C^, equal to zero, which implies that we have a chiral theory 
with either left handed or right handed fermions from the very beginning, we can 
interpret the corresponding energy level crossing of the Dirac field as creation or 
annihilation of fermions. We would like to emphasize that this process takes place 
in real time, in contrast to instanton like effects related to barrier penetration in 
Euchdean space time. 

Finally, we shortly comment on the case of gauge group SU{A). Under as- 
sumption (pSD, we can completely solve the field equations. We get two spinor 
fields with opposite chirality and equal density (C„ = C^), see equations (|83| ) and 



(H), resp. (HD and (||). 

Comparing them to (|8^) we find that the Yang-Mills potential takes values 
only in an su{2) Q) su{2) subalgebra of 5^(4) and, therefore, splits into two parts. 
Each part is coupled to one spinor field of a definite chirality. Thus, in some 
sense, we simply have a doubling of the solution for SU{2). 

As it was mentioned above, equations ( [75| ) and ( |5DD describe the exchange 
of energy between spinor and gauge field. Therefore, in the case of gauge group 
SU{4) we also have energy- level crossing in the evolution of the spinor field. But, 
in contrast to the SU{2)-case, we have - due to constraints ( |69D and (^) - only 
solutions with equal density of left and right handed fermions, i.e. we have no 
chiral solutions. 
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